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Abstract 

The perturbative theory of the nucleation kinetics is analyzed. A 
new improvement is suggested and compared with numerical calcula- 
tions. 

1 Introduction 

The global kinetics of the nucleation process was the subject of interest dur- 
ing the last few decades. Various approaches to the description of nucleation 
kinetics were elaborated in different external conditions pQ - [S] for differ- 
ent systems. One of the most popular is nowadays the perturbation theory 
approach presented in [4J. It was formulated for the external conditions of 
the "decay" type: at the initial moment of time the initial supersaturation is 
created in the system by some external action and later no further external 
influence can be observed - the system evolution occurs only due to internal 
processes of the droplets formation and the vapor consumption by the grow- 
ing droplets. This type of external conditions is rather spread both in the 



experimental investigations and in the theoretical descriptions. The reason 
is evident: on one hand the amplitude value of supersaturation is regulated 
by an external influence, on the other hand the form of the back side of the 
droplets size spectrum and characteristic duration of the nucleation period 
is governed by internal process of droplets formation and growth which gives 
information about the nucleation origin. 

Here we shall analyze the perturbation theory [4] in kinetics of the metastable 
phase decay. The structure of our analysis will be the following: 

• At first we shall see that the recipe given by the perturbation theory 
leads to the non-uniform decompositions 

• Then it will be shown that the perturbation approach in the first ap- 
proximation (it is rather difficult to speak about the further approx- 
imations because even the first approximation can not be calculated 
analytically up to the very end) can be treated as the monodisperse 
approximation 

• A precise solution of evolution equations will be presented and it will 
be compared with the perturbation theory approach. It will be shown 
that the error of the perturbation theory approach will be essential 

• A new version of monodisperse approximation has been proposed. This 
approximation is much more accurate than that given by the pertur- 
bation theory 
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2 Reduction of the balance equation 

The value of supersaturation ( is defined as the ratio 

< = --i 



'<oo 



where n is the number density of molecules in mother phase and n M is the 
number density of molecules in the saturated mother phase. Initial value of 
supersaturation is marked as Co- 

The balance equation in [I] can be written as 

Co 



poo 

-1=A p k g(p,t)dp 
Jo 



at) Jo 

for the distribution function g(p,t) which satisfies the continuity equation 

dg(p,t) _ _ M dg(p, t) 
dt ~ V[) dp 

with initial condition 

,(0 t) = ^ 
m > v(t) 

Here p is the "size" of the embryo which grows with velocity 

c 

v = — 
to 

independent on p, A, k and to are some parameters. Then the number of 
molecules in the embryo of the "size" p will be proportional to p k which 
explains the sense of the first equation as the balance one. 
One can easily reduce this system of equations to 

Co _ x _ A f i u ^nat'))^, 



at) Jo k j v (t> 

with slightly another value of parameter A. 



The case k = is extracted by the possibility of the analytical solution 
of equation ([!]) which can be reduced to the first order differential equation 

with an evident integration. So, we shall discuss this case later as well as the 
case of small fc<Cl which can be solved on the base of the solution at k — 0. 



3 Singular terms in decompositions 

Having extracted the small parameter 

where H is the height of activation barrier in units ksT, one can easily see 
that it approximately equals to the inverse number of molecules in a critical 
cluster. This parameter will be the small parameter of the perturbation 
theory. 

Certainly, after the calculation of all terms in decompositions of the per- 
turbation theory approach and their summation one can get something ac- 
curate. The problem is to get concrete results in frames of approximations 
which can be analytically calculated at least in main features. But already 
in the first approximation of the perturbation theory there appeared an aux- 
iliary function (p^ given by equation 

— = exp(-x ipk) 

MO) = 
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which does not allow analytical solution. So, already the first approximation 
(and all other ones) can not be calculated analytically 

From the last remark it follows that it is necessary to have in the first 
approximation already a good approximation for the real solution. Never- 
theless it will be shown here that the first approximation is not accurate 
enough. Although it was announced in [I] that the decomposition goes on 
the small parameter similar to the inverse number of molecules in droplet, 
the first approximation is far from the real solution. This occurs due to the 
non-uniform character of decompositions. Namely in the balance equation 
(12) in [1] one can see that p k transforms into 

+ + ...)* (3) 

with dimensionless variables 

Cot I h - p 

x = 

e 

playing the role of shifted size and 

t 

T = €Wi 

to 

playing the role of time. All other parameters can be found in [1] . It is evident 
that the r.h.s. is irregular when r goes to zero. From the first point of view 
the limit r —> corresponds to the negligible part of formation of droplets. 
But the careful analysis (see the next section) shows that the nucleation 
period duration has the relative smallness less than e being compared with 
the imaginary time of consumption of the main part of surplus substance. 
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4 Relative smallness of the nucleation period 
duration 



For all k we see that 

Co d(_ r< + + ,,._,um). u , 



( 2 (t) dt Jo" ' v(t') 

for k > and 



Ak [ (t-t'f-^-^P-dt' >0 (4) 
7o V ; v f) V ; 



C 2 (t) eft v(t) 

for = 0. So, we see that the supersaturation is a decreasing function of 
time. 

We exclude the case « 1 from consideration here, since the explicit 
solution for k = has been presented. 

Consider at first the times less than some time t p at which the supersat- 
uration C falls to Co(l ~ e )- Then in the leading term equation ([I]) can be 
rewritten as 

1 " 7^ = (A/v(( )) f\t - t') k im))dt' (6) 
so Jo 

For interval [Co(l — e), Co] one can evidently use an approximation 

/(C(t)) = /(Co)exp(6- 1 ^^) (7) 

so 

Then 

1-^1 = (AT(Co)/«(Co)) A* " t') k exp(- e - 1 (l - (8) 
Co Jo Co 



or 



m = (6- 1 ^J(Co)/^(Co)) A* - t') fc exp(-^(t'))^ (9) 

i/ o 
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for the function 



m = e (i - ^) 

so 



Certainly the last equation can be solved since the renormalization 

t - ( e -M/(Co)/t;(Co)) 1/(/c+1) t 

and 

t > _> ( e -i,4/(c )MCo)) 1/(fe+1) t' 

brings the last equation to equation with no parameters 

if>(t)= f (t - t') k exp(-i)(t'))dt' (10) 
Jo 

Then the function ip is the universal function. Condition ip(t p ) = 1 makes 
the time t p the universal constant. 

For t >t p solutions of ([1]) and ffTUl) practically coincides in their functional 
form. Consider the following equation 



C(t) Jo v(t') 

which is truncated equation (JTJ). Certainly the solution ( of equation 
practically coincides with the solution ( t r of equation ffTTl) for t > t p and is 
less than ( tr for t > t p . So, we get the estimate from above for (. Since 
> I{(2) for two arbitrary £i > ^, we see that we know the estimate 
from above also for the rate of nucleation /. 

One can see that for t > t p the solution of fllip is simple, for integer k 
the integral transforms into polynomial. Then one can easily show that for 
all k except k 1 the solution ffTTj) ensures the rapid decrease of ( tr and 
I((tr)- So, we see the rapid decrease of I(((t)) based on the real solution of 
([T]). Namely, it means that 



• The nucleation period - the period of intensive formation of embryos of 
a new phase is well defined. The lower boundary of this period exists, 
the long tail of the embryos size spectrum is absent. 

• The characteristic relative variation of ( during the nucleation period 
has the order e. 

• The characteristic duration of the nucleation period has the order t p and 
in comparison with the time of essential consumption of the metastable 
phase tfi n it has the smallness e. 

The last property shows that the singularity in is really important 
and can not be neglected. It leads to the very approximate results in the 
first approximation of the perturbation theory. 

One can see here another way to calculate the characteristics of nucle- 
ation. Really, due to the mentioned properties one can spread equation (FlUi) 
to the whole interval [0, oo] of time, get the universal solution and calculate 
universal constants 



Only qi presents all information important for evolution after the nucleation 
period before the coalescense. For integer k the number of qi is finite, for 
arbitrary k one can fulfill decomposition 



The given number of terms ensures practically ideal result, only first three 
qi are necessary. 

This method which is very effective was originally proposed in [5]. 
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5 The monodisperse essence of the first ap- 
proximation 



Now we shall analyze the first approximation in order to see the analytical 
structure. We can consider equation ffTUl) having rewritten it in a form 

M) = (k + 1) f (t - t') k exp{-tfj{t'))dt' (12) 
Jo 

which is more convenient to get some interesting properties. 

The key equation to see the analytical structure is eq. (15) from [I]. One 
can see that the factor p k disappears from the subintegral function and this 
corresponds to the monodisperse approximation of the number of molecules 
in the droplets. Really, in eq. (15) from [I] only the distribution function g 
without p k is integrated. The formal explanation is simple: since the relative 
duration of nucleation period is small one can consider that all droplets are 
formed simultaneously. But this does not work at the nucleation period 
which is the most important because here all droplets are formed. 

Since the monodisperce approximation is adopted it is necessary to get 
an equation for the number of droplets N in this approximation. In the 
mentioned renormalization it looks like 

dN 2 , N k\ 

^ = - 2 -p(-^) (13) 

where 

K k + l J 

It is important that all already formed droplets (and even very small ones) 
have equal contribution in the monodisperce approximation. This leads to an 
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error. The possible way of correction is to use monodisperse approximation 
from pj. It is constructed in the following manner: only droplets formed 
until some moment of time are included in the monodisperce approximation. 
It leads to equation of the following type for the rate of nucleation which is 
proportional to the derivative of N on x: 

I(x)=I(( )exp(-N pr (j)x k ) (14) 

Np r (x) = / I pr (x')dx' 
JO 

Here I pr is the precise solution of ffl2|) . Parameter I means that only the part 
(namely 1/7) of the whole interval produces droplets which are taken into 
account. The value of I can be treated as the characteristic halfwidth of the 
subintegral function p k g, but later more simple definition will be given. 

One can substitute in the last equation the precise solution I pr by the 
value calculated on the base of the same monodisperce approximation. Then 
we come to the self-consistent monodisperce approximation 

I S c(x)=I((o)exp(-N sc (j)x k ) (15) 

N sc (x) = f I sc {x')dx' 
Jo 

It is possible to continue the simplification of the monodisperse approxi- 
mation. One can see that at the nucleation period the argument of exponent 
in expression for I increases very rapidly. So, it is quite possible to assume 
that for x at nucleation period the argument x/l corresponds to practically 
ideal situation when the action of droplets is negligible even on the rate of 
nucleation. It can be done when I 1. Then we come to 

ry ry> ry> 

tf(j)«*i(yW(Co)j 
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We see that the functional form of / here is 



k+l 



I{ x ) = J (Co) exp( 



(16) 



and it coincides with the functional form in the iteration solution of f|T2|) . 
The iteration procedure is defined in [3] as 



It is known that this is rather good approximation for the rate of nucleation 
aat least at k = 3/2, 2, 3 typical for nucleation. 

The full coincidence will be if we put I = k+l. This is the most reasonable 
choice of I. One can see that here / attains big values as it has been assumed. 

Certainly one can choose / in a more sophisticated style. The problem 
is to write and to solve some algebraic equation on I which can be done in 
more precise approximations. It is not a real problem and here it is quite 
sufficient to use the simplest I = k+l. In any case it can be refined. 

Now it is worth to return to the problem of analytical calculation of the 
mentioned approximations. Now only (ITS!) can be calculated in analytical 
way. The way to calculate other approximations is the following: 

• Consider ( TLB"]) as the base approximation. Having presented / as I = 
Iij- we decompose exp(ln/ — ln/i) into series of argument (In I — In Ji). 
It is worth doing because both / and I\ have exponential form. Then 
we come to expression of a type 

rx 

/ exp(— x k+1 ) * Polinomial(x)dx 








(17) 



and leads to 




.fc+i 
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which can be easily reduced to the sum of ^-functions or f£ exp(—x l )dx. 
The problem is solved. 

The same method can be applied to calculate further iterations in the 
iteration procedure (IT7|) . 

It seems that new approximations will work better than the old ones. 
Now we shall test all of them. 

6 Discussion 

Here we shall investigate in details all cases which were calculated in [I], 
namely k = 0, 1/2, 1,2. We add here the case k = 3 because it corresponds to 
formation of embryos under the free-molecular regime in a three-dimensional 
space which is the most natural case. 

The case k = has the analytical solution as the first order differential 
equation. It is shown in figure 1. One can see that here the rate of nucleation 
as a function of time (or of coordinate x) is shown. 
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Figure 1 
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The rate of nucleation as a function of x forms the droplets size spectrum. 
One can see that here it is rather long. The tail plays the main role, the 
number of droplets is infinite at least under the exponential approximation 
of the nucleation rate. 

So, here appears a problem to calculate the total number of droplets. The 
exponential approximation fails in this question and this opens a problem. 
Fortunately, this problem has a very simple solution. 

One can see that here it is extremely simple to calculate the main charac- 
teristic of nucleation - the total number of droplets. Really, the total number 
of droplets N tot is simply the amount of surplus substance Co^oo divided by 
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the final number of molecules in the droplet Vf in : 



^ 'fin 



Here we see the error of the last model because for k = the rate of growth 
is zero and the droplets can not become supercritical ones. So, the necessity 
of external parameter u Hm is evident. 

The case of small k is the most dangerous because on one hand there is no 
such balance relation as the previous one and on the other hand all droplets 
play practically equal role in vapor consumption and the number of droplets 
appeared at the fallen supersaturation is essential. 

Consider now the power k = 1/2. The rates of nucleation in different 
approximations in this case are drawn in figure 2. 



Figure 2 
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The rate of nucleation is drawn. Here one can see two point curves and 
two line curves. The upper point curve is the precise solution. Solution with 
decomposition of kernel can not be separated from the real solution in scale 
of a picture - the error is negligible. The upper line curve is the approx- 
imation of the perturbation theory. It is rather close to the real solution, 
but still the error exists. The lower line curve is thick - two monodisperse 
approximations IHM and (fl5|) are drawn here. They can be hardly separated. 
The approximation of the first iteration (1161) is the lower point line. One can 
see that here the error of the last approximation is essential. 

Consider now the case k — 1. This case can be solved analytically, nu- 
merical results are drawn in figure 3. 
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Figure 3 
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The rate of nucleation is drawn. Here one can see four curves. The upper 
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point curve is precise solution. The upper line curve is approximation of the 
perturbation theory. The lower line curve is thick - there are two monodis- 
perse approximations ffl4|) and ffl5l) which again can not be separated. The 
lower point curve is approximation of the first iteration. All approximations 
work well but the perturbation theory approximation has the maximum er- 
ror. Later this tendency will be stronger. 

The next case is k — 2. It is shown in figure 4. 
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Figure 4 
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The rate of nucleation is drawn. Only two curves can be seen here. The 
curves corresponding to the precise solution and approximations (EE]) , (|T5l) 
and ffl6|) is one thick line with points on it. Another line is the approximation 
of perturbation theory. Although it has an evident error it seems that with 
k becoming great this error will be small. But it is no more than illusion. 
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For k = 3 we have the following numerical results which are drawn in 
figure 5. 
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Figure 5 
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The rate of nucleation is drawn. Again one can see two curves here. The 
curve without points corresponds to the approximation of the perturbation 
theory. The curve with points corresponds to precise solution and all other 
approximations. One can see that being compared with the case k = 2 the 
approximation of perturbation theory goes away from precise solution. This 
tendency takes place because of a wrong number of essential droplets in the 
actually monodisperse approximation of the perturbation theory. 

Now one can calculate the errors of new approximations. Since the errors 
of the perturbation theory can be seen directly in pictures there is no need 
to discuss it. The errors of approximations for k — 1, k — 2, k — 3 are drawn 
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in figure 6. 




One can see here three pairs of point and line curves. The line curves 
are approximations f|T6|) for different k, the point curves are approximations 
(TH|) , (fT5j) (they practically coincides for all k with the reasonable accuracy of 
calculations). The lowest point and line curves coincides (this corresponds to 
k = 3) and it means that for k > 3 there is no advantage of approximations 
(|T4|) instead of (fT6l) . The upper curves correspond to = 1, the middle 
curves - to k = 2. One can see that for k = 1 the advantage of (JHJ), (|T5l) is 
essential. 

The approximation f|T6|) can lead to some important consequences. Hav- 
ing drawn in figure 7 the approximation f|T6l) for all considered = 1/2, 1, 2, 3 
we see that there is a focus - a point where all curves get together. 
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Figure 7 
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This fact can be proven analytically by the simple differentiation. Really, 
the form of the first iteration is extremely simple ~ exp(x ( - fc+1 - ) ) and at x — 1 
the rate of nucleation falls e-times. 

Since the precise solution is not far from the first iteration, one can see 
the approximate property. It is shown in figure 8 where the precise solutions 
for k — 1/2, 1, 2, 3 are drawn. 
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Figure 8 
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This allows to speak about the property of the precise common length of 
nucleation under the external conditions of the decay type. 
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